Intermittency and non-Gaussian fluctuations 
of the global energy transfer in fully developed turbulence 



CO 

o 
o 

CN 



CO 

o 



^3 

C 

o 

CN 
> 
CO 

o 
in 

CN 



i 

C 

o 
o 



X 



B. Portelli, P.C.W. Holdsworth and J.-F. Pinton 

Laboratoire de Physique, Ecole Normale Superieure, 
46 Allee d'ltalie, F-69364 Lyon cedex 07, France 

We address the experimentally observed non-Gaussian fluctuations for the energy injected into 
a closed turbulent flow at fixed Reynolds number. We propose that the power fluctuations mirror 
the internal kinetic energy fluctuations. Using a stochastic cascade model, we construct the excess 
kinetic energy as the sum over the energy transfers at different levels of the cascade. We find 
an asymmetric distribution that strongly resembles the experimental data. The asymmetry is an 
explicit consequence of intermittency and the global measure is dominated by small scale events 
correlated over the entire system. Our calculation is consistent with the statistical analogy recently 
made between a confined turbulent flow and a critical system of finite size. 
PACS : 47.27.-i,64.60.Cn 



Intermittency is a well established feature of turbulent 
flows. It has been observed in a large number experi- 
ments and incorporated in many models In this pa- 
per we attempt to link the statistics of the fluctuations of 
the global energy consumption Pj in a turbulent flow to 
the intermittency of the energy transfer. Our motivation 
stems from the empirical observation [2j that the proba- 
bility density function (PDF) of power consumption Pj 
for a confined turbulent flow driven at constant Reynolds 
number Re has a non-Gaussian form, asymmetric with 
an exponential tail 0, 0] , extremely similar to that ob- 
served for the order parameter fluctuations in a finite size 
equilibrium system at criticality. Very similar distribu- 
tions can be observed for global quantities in other ex- 
perimental [E IE IE El an d model correlated systems 0] • 
These observations have generated recent interest in di- 
verse aspects of global fluctuations, fo r sy stems both in 
and out of equilibrium [III 111 EE Q El 13 We have 
previously argued that the observation of similar fluctu- 
ations in the global quantity for these radically different 
correlated systems implies an analogy at the statistical 
level 0. Specifically, we have shown 0] using the low 
temperature phase of the 2D-XY model, that a key in- 
gredient for the occurrence of such non-Gaussian fluctua- 
tions is the action of modes with a diverging distribution 
of spin wave stiffness across spatial scales. In order to ap- 
ply these findings to the statistics of power injection in 
a confined turbulent flow, we consider turbulence in the 
framework of a stochastic cascade across scales. The flow 
state, at a statistical level is viewed as gas of indepen- 
dent fluctuations at all scales I between the injection and 
dissipative lengths, rj and L. The variance of the fluctu- 
ations of these modes of energy transfer diverge at small 
scale, which is the well accepted feature of turbulence, 
called intermittency . The divergence across the scales 
is analogous to that observed in the model critical system. 
In this paper we therefore show, using a phenomenolog- 
ical model, that the observed non-Gaussian fluctuations 
of the global power injection are a natural consequence 
of intermittency. 



Experimental data, reproduced from reference 5], are 
shown in figure 1. The flow is driven by two concentric, 
counter rotating discs, with cylindrical geometry, rotat- 
ing at constant frequency f2. The probability distribution 
for the fluctuations in the power, Pi{t), needed to drive 
the flow, is characterized by a negative skewness, with 
an approximately exponential tail for fluctuations below 
the mean and a much more rapid fall-off for fluctuations 
above the mean. The energy balance is given by U 



dE(t) 
dt 



Pi(t)-P D {t) 



(1) 



where E(t) and Pd(*) are respectively the instantaneous 
kinetic energy and dissipation by viscous forces 0, of 
the confined fluid. A time average of equation (1) yields 
the condition of stationarity: Pj = Pp = Mcq, where 
we define eo, the power consumption of the flow per unit 
mass, and M, the total mass of the fluid. 

In the absence of a microscopic description of the sta- 
tistical properties of the flow, a phcnomcnological model 
is required to explain the observed fluctuations about this 
mean value: firstly, we propose that the fluctuations in 
Pi(t) are related to fluctuations in the internal kinetic 
energy of the flow 



P/(t) = Me - £{t) 



(2) 



Here, £ it) — (E(i) — E)/t is the difference between the 
instantaneous and mean kinetic energy, normalized by t, 
an integral time scale characteristic of the energy injec- 
tion and of the response of the driving mechanism (in- 
ertia). That is, we propose that the torque required to 
shear a fluid at constant rate Q is a decreasing function 
of the internal kinetic energy of the fluid. We predict, 
therefore, that the constraint, 17 = const., forces the fluc- 
tuations in Pi(t) to be a mirror reflection of the fluctu- 
ations of the internal energy of the flow. If constrained 
with constant torque one would consequently expect re- 
versed power fluctuations with positive skewness, as also 
observed experimentally We do not suppose that 
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mechanical structures drive the disks, simply that the en- 
gagement of the disks is reduced with increasing kinetic 
energy of the underlying flow. An excess of kinetic energy 
should be repartitioned over the entire flow and felt si- 
multaneously at the two disks, which would then decrease 
their engagement in the fluid, independently of their ro- 
tation direction. This scenario provides a stabilizing loop 
through which the turbulent steady state is maintained 
and it is consistent with experimental observation that 
the time variation of the power injected into the upper 
and lower discs displays positive correlations: power fluc- 
tuations of all amplitudes are felt quasi-simultaneously at 
both disks, despite their rotation in opposite directions 
(See Fig. la of ref. @|). 

Secondly, we build £(t) as the excess of transferred 
energy at time t, summed over contributions from all 
scales in the flow. For this we introduce nt(r,t), the 
instantaneous real-space energy flux: 

Mr,t) = -~V*- (\6eu(r)\ 2 8 e u(r)) , (3) 

which we average over the flow volume to obtain the en- 
ergy flux at scale £ 

= ~ J d 3 r n t (f,t) , (4) 

these variables could be computed in a direct numeri- 
cal simulation. Our phenomenology concerns the fluc- 
tuations away from the steady state at all scales £: 
{^i{t) — en). We adopt a statistical point of view rather 
than a dynamical one and we make the assumption that 
the (iif — en) are statistically independent stochastic vari- 
ables 19] (one should think of this as a non-linear change 
of variables rather than a physical tree structure replac- 
ing the real flow) . We then define the excess energy £ as 
the sum over contributions from all scales: 

£ = M^fa - eo) . (5) 

i 

This is a strong approximation, whose justification rests 
with the results it produces, which are in striking agree- 
ment with experiment-see Fig. 1. Physically, our model 
implies that the energy injection at a given time takes 
into account the entire structure of the energy flow at 
that time. When 7i^(t) > en, energy is built up at scale 
I and the overall energy excess £(t) is the sum of such 
effects over the entire range of scales. 

Equation (JSJ, together with our assumptions allows 
one to construct the PDF of £ as the convolution of prob- 
ability densities for the elements (ni — en): 

n(e) = Y[p e (n-eo) , (6) 

where pt{iri) is the PDF for the fluctuations of tt£. Once 
the pe(ire) are chosen, the above expression allows for the 



numerical computation of the distribution for the global 
quantity. In fact, one can express the PDF for the nor- 
malized variable (f> = £ ^ £> as 

r+°° f ]k 

m = J ^e ik **(k) , (7) 

where the characteristic function 'i'(fc) is the product of 
the characteristic functions of each of the microscopic 
distributions and the net variance a 2 is the sum of the 
microscopic variances. Then, using equation (2), the 
PDF for the fluctuations of the normalized power input 

(0 = fi=E)fen(0)=ii(-0). 

The phenomenology proposed here cannot be devel- 
oped within the context of Kolmogorov's 1941 (K41) the- 
ory of turbulence 0, as it ignores fluctuations around 
mean values of energy injection and transfer. We there- 
fore develop equation (0 in the context of Kolmogorov's 
Refined Similarity Hypothesis 0, E3- ^ allows for 
anomalous scaling of the energy transfer at scale £: 

(7r?)«cg (e/L) TM . (8) 

The spectrum r(q) is related to the experimentally ob- 
served non-linearity in the scaling exponents of the 
longitudinal velocity increments: (Siu q ) oc £^^ cx 
(Tr q / 3 ) (£/L) q/3 so that ((q) = q/3 + r(q/3). 

In the following we use two different forms for the 
microscopic distributions, compatible with the experi- 
mentally determined exponents for the velocity struc- 
ture functions £(<?). We show that intermittency is nec- 
essary to obtain asymmetric non-Gaussian distributions 
with negative skewness. However, let us first return to 
some details of the model: the {L/rf) 3 correlated degrees 
of freedom implicated in the flow Q are transformed 
into an ensemble of N scales. At each level the system 
is described with a resolution length scale £ n and the 
scales are separated by a ratio A = (£ n -i/£ n ) 3 such that 
X N = {L/rf) 3 . The Reynolds number and the number of 
levels are related through the definition Re = (L/rj) 4 / 3 = 
\ 4N / 9 , Note that the variable A is a free parameter in 
the theory. In the original Kolmogorov-Obukhov theory 
(K062) of intermittency 0, the energy cascade, being 
a multiplicative process with a large number of steps, 
is assumed to have log-normal statistics. In this case 
p n (log TT£ n ) is a Gaussian distribution whose variance can 
be that of K062 or that proposed by Castaing et al. [2l| . 
As one assumes scale invariance in K062, the variance 
is <j2 = e 2 [(£ n /L) T ( 2 ) - l] . This distribution leads to 
the quadratic spectrum ((q) = q/3 + pq(3 — g) / 18. The 
value of the intermittency parameter (i = — t( 2) that fits 
best the experimental data is p — 0.21 Another 
possibility, which gives a very good fit to the velocity 
intermittency exponents is a \ 2 distribution: 

Pn(ire n ) = N e n ' ir ? n ~ ex P ( _a eii n Tre n )- (9) 
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In this case, we require that the variance of p n (%£ n ) be 
that of K062, in addition to the conditions of normal- 
ization and constant mean. The three conditions deter- 
mine N# n , ai n , vi n . From the expression for the moments 
(< >= £qF(<7 + Vi^)jv\ T{vg n )) it is straightforward 
to show that the spectrum r{q) is that of K062 in the 
limit l n — > L. Even if x 2 statistics does not strictly 
allow for scale invariance for r(q) and £(q) in the inertial 
range, the corrections are very small. The main advan- 
tage of the x 2 distribution is that it allows a straight- 
forward calculation of the convolution product in equa- 
tion 0: 

N 

log*(fc) =ike /a--^ f v in log(l + ik/(aa tn )) (10) 

n=l 

— . ' ' ■ 



c 




FIG. 1: PDF of global energy transfer for experimental data, 
taken from Pinton et al. [f| (symbols) and from cascade mod- 
els of intermittency. The cascade models are: x 2 -model, 



Re 
Re 



- 10 5 and A = 2.0 (solid line) and lognormal model, 
10 s and A = 1.08 (dashed line). 



In figure 1 we show the experimental data published in 
ref || together with the PDF for n(6*), calculated from 
equation (5), using both log- normal and x 2 models. They 
both capture the essential features of the experimental 
data: (i) the fluctuations in 9 are strongly non-Gaussian, 
despite the high value of the Reynolds number, i.e. of the 
large number of contributions in the sum JSJ; (ii) the dis- 
tribution is skewed towards negative values; (iii) the tails 
of the PDF towards negative values are almost exponen- 
tial, although this seems better verified for the \ 2 micro- 
scopic distribution than for the log-normal. This behav- 
ior is extremely robust against variations of the Reynolds 
number, at fixed A. In this case, increasing the Reynolds 
number leads to an increase of the number of cascade 
steps TV (as logi?e) and, at the same time, to an increase 
of the asymmetry in the shape of the microscopic distri- 
bution at the smallest scales. In figure 2 we show the 
evolution of the distribution with Re, for A = 2.0, for the 
X 2 microscopic distributions. The curves vary only slowly 
with Re. The asymmetry becomes more pronounced as 
Re increases, as can be seen from the skewness, 7, which 



varies from —0.64 to —0.97 for Re varying from 10 4 to 
10 7 . Hence, as in the experiment, there is no evidence 
that 11(6*) will reduce to a Gaussian in the limit of infinite 
Reynolds numbers, and therefore infinite N. Given the 
uncertainty in the experimental results and the crudeness 
of our model, this qualitative agreement seems extremely 
encouraging. 

As A is a free parameter in the problem, it has been 
chosen in figure 1 to give a good fit of the experimental 
curves for each class of microscopic distribution. For a 
fixed value of Re, varying A does not change the shape 
of the global PDF. Using the x 2 distribution, we observe 
a slight increase of the skewness, from —0.65 to —0.83 
when A is increased from 1.5 to 2 (for Re=10 6 ), a typical 
range of values used in cascade models . 




FIG. 2: PDF of global transferred energy for Re varying be- 
tween 10 4 and 10 7 (legend). The microscopic distributions 
are \ 2 ones and A = 2.0. 



However, there seems no a priori reason to fix A and 
while we do not rule out a weak variation, we can fix a 
universal form for the PDF, as proposed in Q by letting 
A vary with Re. Support for this proposition comes from 
the evolution of the ratio cr/Pj with Re. Experimentally, 
one observes a very slow decrease that may be fitted by a 
power law dependence with exponent a ~ 0.3 In our 
computations at fixed A, we find the ratio to be slightly 
increasing, in disagreement with experiment. If, instead 
of fixing A we impose universality, we find a decreasing 
ratio, with a weak power law dependence with exponent 
a ~ 0.1. 

Here, the small scales play a major role due to the dis- 
persion in the amplitudes and as in the finite size critical 
system j^, one might expect the global PDF to retain 
some characteristics of the microscopic measure, with the 
experimentally observed asymmetry reflecting the form 
of Pn (""£„)• This is indeed what we observe and the sym- 
metry of the global PDF is only restored if the individual 
elements are made to have equivalent weight. This is il- 
lustrated in figure 3, where we show a series of PDFs 
produced at fixed Re and A, with r(2) varying between 
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the experimental value, 0.21 and zero. As r(2) is reduced 
the PDF becomes more symmetric and, for r(2) = it is 
well represented by a Gaussian. This point is quite im- 
portant: intermittency is absolutely necessary to obtain 
an asymmetric and non-Gaussian PDF for the global en- 
ergy transfer. However this conclusion only holds for the 
PDF of energy transfer. In the case of pressure fluctua- 
tions Holzer et al. |24| have shown that intermittency is 
irrelevant since for the pressure one integrates Siu 2 which 
gives a dispersion of variance < Siu 2 >~ (eal) 2 ^ 3 in the 
strict KA1 description. We remark that, with the phe- 
nomenology proposed here, negative skewness is a direct 
consequence of kinetic energy fluctuations with positive 
skewness, as are observed in microscopic models with dis- 
sipation [HE1E3. 

The importance of the small scale can be gauged by 
recalculating Z , excluding levels from either the top or 
the bottom of the cascade. We find that with removal of 
cascade levels near the dissipative scale the distribution 
rapidly approaches a Gaussian form, while removal of lev- 
els from the top of the cascade leaves the PDF essentially 
unchanged. The entire global measure is therefore seen 
to be influenced by a relatively small number of statis- 
tically independent cascade levels describing the system 
at small scales. If a universal PDF for the global mea- 
sure is indeed an experimental reality, then the disper- 
sion in amplitudes must diverge with Reynolds number, 
in such a way that the universal distribution is a "limit" 
function, valid even as Re — * oo. These arguments lead 
us to a striking physical conclusion concerning the na- 
ture of a turbulent flow: given the huge number of small 
scale objects present in the flow ((i/77) 3 = Re 9 / 4 ~ 10 9 
for Re = 10 4 ); if the PDF for the global measure is 
non-Gaussian, it can only mean that these objects are 
strongly correlated in space and time. This point seems 
consistent with recent measurements of long time tem- 
poral correlations in the dynamics of Lagrangian tracer 
particles in a fully turbulent flow |25| . 

We finally address the point that the PDF of the in- 
jected power has a Gaussian shape for open (non con- 
fined) flows Q. The distribution of the global quantity 
thus depends on the overall size of the flow, although it 
is known that small scale intermittency characteristics 
do not; for instance, the scaling exponents of the longi- 
tudinal velocity increments are identical in open or con- 
fined flows |22(. We propose that in non-confined flows 
several uncorrelated cascades occur simultaneously. In 
this case, the probability for the overall energy transfer 
is given as the convolution of J\f functional forms of the 
type H(0) computed above. This tends rapidly to a Gaus- 
sian. Experimentally, this could be tested by changing 
continuously the ratio of the disk to cylinder radius in a 
'washing-machine' experimental setup. 

The experimental data shown in figure 1, published 
in reference have been obtained with R. Labbe. We 
gratefully acknowledge many fruitful discussions with M. 








FIG. 3: Variation of the global PDF with the magnitude of 
the intermittency parameter |r(2)| for Re — 10 6 and A = 2.0. 
The microscopic distributions are \ 2 statistics. 
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